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Abstract 

We give a detailed analytical description of the global dynamics of a point mass moving on a 
sphere under the action of a logarithmic potential. After performing a McGehee-type blow-up 
in order to cope with the singularity of the potential, we investigate the rest-points of the flow, 
f~* > ) ■ the invariant (stable and unstable) manifolds and we give a complete dynamical description of 

£SJ ' the motion. 
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Introduction 

Q. 

Topologically, two dimensional Riemann surfaces with constant (Gaussian) curvature K are clas- 
sified into three categories: Euclidean spheres, S 2 (K > 0); Euclidean planes, E 2 {K = 0); and 
hyperbolic planes H 2 (K < 0). Among them, S 2 and E 2 are more familiar and come out very 
often in practice. For example, the mechanics of thin fluid layers on S 2 provides a global model of 
a planetary atmosphere, and on E 2 its local approximation. 

In this paper we analyze the motion of a point particle moving on a sphere under the action of 
a logarithmic potential. Two are the main reasons for the choice of this particular potential. First, 
it arises in different physical scenarios: such as in models of astrodynamics, [ToTr97] , |StFo03j : in 
£NJ , the dynamics of a charged particle in a cylindrically symmetric electric field [Hoo63 and in the 

mathematical theory of vortex filaments of an ideal fluid NcwOl , [KcPoVc03 . The second reason 
relies on the fact that the logarithmic potential V(x) — — log(|a;|) could be considered as a limit case 
for a — > of the homogeneous potentials V a — \x\~ a and, while the latter have been extensively 
studied by different authors, the former has not been so deeply investigated. In particular one 
could be interested to know if (and how) some features regarding for instance the regularization of 
collisions, the minimality properties of the solutions, the stability character, may be extended from 
the homogeneous to the logarithmic potential case. Results in this direction have been achieved for 
instance in [Cas09] . |BaFeTe08j . |ToTr97j . 

In addition, we consider a sphere, rather than the classical two or three-dimensional Euclidean 
space, as the configuration space. Our goal is to understand which aspects of the dynamics are 
affected if the geometry of the underlying space changes, or equally well, what survives of the planar 
dynamics if one considers a curved manifold. 

From a dynamical point of view the most interesting feature, and the hardest obstacle for a 
full understanding of the motion, is played by the presence of the singularity in the potential 
function. Indeed, as it often happens in celestial mechanics, the singularities are the source of a 
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complicated dynamics and sometimes they are even responsible for a sort of chaotic motion. From 
the mathematical point of view the singularities represent a severe technical hurdle to overcome and 
different techniques have been proposed to regularize the vector field, , mainly for the homogeneous 
potential case |LeC20j . |McG74j . |Eas71| . |BeFuGr03], |Cas09j and [CaTe] . 

This paper is inspired by the recent a work IStFo03] that studies the planar motion of a point 
mass subject to a logarithmic potential in an astrodynamic context. To overcome the singularity 
of the vector field we adapt to our problem the celebrated McGehee transformation, a regularizing 
change of variables currently popular in the field of Celestial Mechanics and first introduced in 1974 
by McGehee [McG74] to solve the collisions in the collinear three-body problem. 

The McGehee transformations consist of a polar type change of coordinates in the configuration 
space, together with a suitable rescaling of the momentum. In this way the total collision is blown- 
up into an invariant manifold called total collision manifold over which the flow extends smoothly. 
Furthermore, each hypersurface of constant energy has this manifold as a boundary. By rescaling 
time in a suitable way, it is possible to study qualitative properties of the solutions close to total 
collision, obtaining a precise characterization of the singular solutions. 

The McGehee transformation are usually applied to the case of homogeneous potentials but, as 
shown in [StFo03 and as it will be manifest throughout this paper, with slight modifications they 
give interesting results even in the presence of a logarithmic potential. In fact, although the lack 
of homogeneity of the logarithmic nonlinearity breaks down some nice and useful properties of the 
transformation, it is still possible to regularize the vector field and therefore it is still possible to 
carry out a detailed analytical description of the rest points, of the invariant manifolds, and of the 
heteroclinics on the collision manifold. 

The paper is organized as follows: first we introduce some basic notions about the Hamiltonian 
formulation of the co-geodesics flow on a general Riemanian manifold, then in Section [2] we restrict 
to the case of the sphere and we formulate the equivalent co-geodesics flow through the stereographic 
projection. In Section [3] we introduce the singular logarithmic potential and we write the equation 
of motion. Section 2] deals with the in-deep study of the dynamical system: we regularize the 
singularity of the potential with the modified McGehee technique, and we provide an analysis of 
the flow on the collision and the zero velocity manifolds. Section [5] concerns the global dynamics 
and we rephrase the results in terms of the original motion on the sphere with untransformed 
coordinates. 
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1 Preliminaries 

Let M be a Riemannian manifold, namely a smooth n-dimensional manifold M endowed with a 
metric given by a positive definite (non-degenerate) symmetric two-form g. We denote by D the 
associated Levi-Civita connection and by ^ the covariant derivative of a vector field along a smooth 
curve 7. Let / be an interval on the real line and let V be a smooth function defined on / x M. 

Definition 1.1. A perturbed geodesic abbreviated as p-geodesic is a smooth curve 7: I — > M 
which satisfies the differential equation 

(1) ^7'(t)+W(i,7(*))=0 
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where \7V denotes the gradient of V(t, — ) with respect to the metric g. 

Remark 1.2. From a dynamical viewpoint, the data (g, V) define a mechanical system on the 
manifold M, with kinetic energy hg(v,v) and time dependent potential energy V. Solutions of the 
differential equation Jl} correspond to trajectories of particles moving on the Riemannian manifold 
in the presence of the potential V. If the potential vanishes we get trajectories of free particles 
and hence geodesies on M. This motivates the suggestive name, "perturbed geodesies" in the case 
W 7^ 0. Moreover, if the potential V is time independent, modulo reparametrization, perturbed 
geodesies become geodesies of the Jacobi metric associated to (g,V): indeed the total energy 

is constant along the any trajectory 7 thus, whenever V is bounded from above, the solutions of (JTJ) 
with energy e > sup m6M V(m) are nothing but reparametrized geodesies for metric [e — V]g on M 
with total energy one fAbMa78^ . 

Denoting by (</,..., q n ) a local system of coordinates on M, equation ([TJ reduces to 

. , ..dV 
q +T ]k q J q =-9 3 q~, 

where, as usual, g^ = and T l - k are the Christoffel symbols. 



Geodesic flow as Hamiltonian flow 

The geodesic flow turns out to be a Hamiltonian flow of a special Hamiltonian vector field defined on 
the cotangent bundle of the manifold. The Hamiltonian depends on the metric on the manifold and 
it is a quadratic form consisting entirely of the kinetic term. The geodesic equation corresponds to 
a second-order nonlinear ordinary differential system. Therefore by suitably defining the momenta 
it can be re- written as first-order Hamiltonian system. 

More explicitly, let us consider a local trivialization chart of the cotangent bundle T*M 



T*M 



U x 



where U is an open subset of the manifold M, and the tangent space is of rank n. Let us denote 
by (qi, q2, . . • , q n ,Pi,P2, ■ ■ ■ ,Pn) the local coordinates on T*M and introduce the Hamiltonian 

(2) H : T*M ->■ K : H(q,p) = ^(x)p lP] . 

The Hamilton- Jacobi equations of the geodesic equation with respect to the metric g can be written 
as 

dH at ^ 
-dp- =9(x)Pj 

dH 1 dg jk 
"dq~^'^W PjPk ' 
The second order geodesic equations are easily obtained by substitution of one into the other. The 
flow determined by these equations is called the co-geodesic flow, while the flow induced by the first 
equation on the tangent bundle is called geodesic flow. Thus, the geodesic lines are the projections 
of integral curves of the geodesic flow onto the manifold M. 

Being the Hamiltonian H time-independent, it is readily seen that the Hamiltonian is constant 
along the geodesies. Thus, the co-geodesic flow splits the cotangent bundle into level sets of constant 
energy 

M E = {(q,p)eT*M: H(q,p) = E}, 
for each energy E > , so that 

T*M = \J M E - 

e>o 
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Now let g, h be two Riemannian metrics on M in the same conformal class; namely there exists a 
positive and smooth function A = X(q) of the coordinates such that 

g ij = X h" 

or equivalently g = A _1 ft.. From the definition (|2|) it follows that a scaled co-geodesic Hamiltonian 
function corresponds to a conformal change of the metric. In fact if H g and Hh denote the Hamil- 
tonian co-geodesic functions and if g, h are in the same conformal class then it immediately follows 
by © that 

Hg(q-p) = XH h (q,p). 

As a consequence, Hamilton's equations with respect to this two Hamiltonian functions are related 
as follows 



(3) 



P 



d p H g 



d q H g 



Xd p Hh + HhdpX — XdpHh 
—Xd q Hh — Hhd q X 



where the last equality in the first equation comes by the fact the the function A only depends on 
<?• 

In the following we consider a perturbed-geodesic flow on a sphere, thus it is worth to write 
down explicitly the free geodesic flow when the manifold M is a surface of revolution in M 3 . 

Denote with (x, y) the Cartesian coordinates of M 2 and consider the function (p : U C K 2 — > M 3 
given by ip(x, y) = (f(y) cos x, f(y) sin x, g{y)) 

U = {(x, y) e M 2 : < x < 2tt, y < y < yi}, 

where / and g are differentiable functions, with f'(y) 2 +g'(y) 2 ^ and f(y) ^ 0. Thus ip(x, y) is an 
immersion and the image f(U) is the surface generated by the rotation of the curve (0, f(y),g(y)) 
around the z axis0 The induced Riemannian metric g = (gij) in the (x, y) coordinates is given by 

ffii =f 312=0 g 22 = {!'f + {g'f- 
From the Hamiltonian function associated to the geodesic flow is given by 



1 ( 1 



H(x,y,p x ,p y ) = - —p. 



1 



2 V/ 



/' 2 + g' 



and the co-geodesics flow reads as 



(4) 



or, equivalently, as 



r- 



i 



y = 

Px 

Py 



f' 2 + g' 2 
o 

//' 2 
^Px 



Py 



f'f' + g'g" 

(/' 2 +5' 2 ) 2 



2//'. . 



f 2 , 
ff 

f' 2 + g' 



ff" 



g'g" 



(y) 2 



2 The stereographic projection of the sphere 

It turns out that transformations of the McGehee type may be devised without too many difficul- 
ties for equations which are written in Cartesian coordinates on a plane. Therefore, rather than 



^^Here we are considering the Euclidean space equipped with Cartesian coordinates whose axis are labeled as x, y, z 
according to the ordering induced by the canonical orthonormal basis of M 3 . 
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Figure 1: Mutual positions of the sphere (R = 1) and the plane. 



attempting to work directly onto the sphere, we felt it would be more easy (and more clear) first to 
project the dynamics on a stereographic plane, and then to remove the singularities of the resulting 
equations. 

We work on a two-dimensional spherical surface S of radius R and center at the point C — 
(0, 0, R), namely 

S := {(re, y, z) e K 3 : x 2 + y 2 + (z - R) 2 = R 2 }. 

where (x, y, z) are the Cartesian coordinates in R 3 (figure [1]). We shall call north pole and south 
pole the point N := (0, . . . , 0, 2R) G S and its antipodal S := (0, . . . , 0, 0) G S, respectively. Note 
that the sphere is tangent at the origin to the plane {z = 0}, that we identify with M 2 . Next we 
introduce the stereographic projection 



tts-S\ {N} 

p 



p, 



defined by requiring that the three points N, P, P are collinear. By a straightforward calculation it 
follows that the map 7rg is given explicitly by 



(5) 



2R 



2R- z 

We use slightly non-standard angular coordinates <f> and 6 for the spherical surface: 

• (j) € [0, 27r) is the usual polar angle of the projection of P onto the plane z = 0; 

• 9 € [0, 7r] is the angle between the segment PC and the negative direction of the z-axis. 

A generic point P = (x, y, z) on the sphere in these coordinates has a local parameterization given 
by 



X 




sin 


9coscj) 




y 


= R 


sin 


9sm<p 




z 




1 - 


cos 8 





Of course the map is a diffeomorphism of class ' 
its is defined as : 

sin 9 cos 4- 
R sin 9 sin <fi 
1 — cos 9 



In these coordinates the stereographic projection 
2R 



1 + cos t 



sin v cos ( 
sin 9 sin t 
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We recall that if M C K 3 is a portion of a regular surface represented in Cartesian local coordinates 
by the vector equation 

P(u, v) := + x(u, v)i + y(u, v) j + z{u, v) k 

then dP = P u du + P v dv for P u — (x u , y u ,z u ) and P v = {x v , y v , z v ) and hence the metric is given 
by ds 2 — \\dP\\ 2 . With the above parametrization of the sphere it follows that 



R (cos 8 cos <f), cos 9 sin (j>, sin t 



P^ = R(— sin 9 sin <f>, sin 9 cos </>, 0). 



Denoting by g and gs respectively the Riemannian metric on the sphere S and the metric on the 
plane induced by the stereographic projection, we have 



g:=R 2 



sin 



gs ■= 



(l + costf) 5 



As a consequence of the above calculation the following result holds. 

Lemma 2.1. The manifolds {S,g) and (M. 2 ,gs) are in the same conformal class, where R 2 denotes 
the Alexandroff compactification ofM. 2 . 



Co-geodesic flows 

Let G and Gs be the matrices corresponding to the inverse of g and gs respectively; thus we have 



„ D _ 2 /sin" 2 9 0\ , _ (l + cos(9) 2 

G = R \ „ „ , and Gs = 



V 



1 



G. 



We denote by S = S \ {N, S}, the sphere minus the north and the south pole, and by T*S its 
cotangent bundle where it is well-defined the Hamiltonian function 



H, 



geod 



t*s -> K 



Here q := (0,6*) are the positions and p := (p$,pe) are the momenta. With this choice the 
Hamiltonian function is given by 



Hgeod(4>: 0',P<t>,Pe) 



1 



1 



2R 2 V sm- 

and, as a particular case of d3|, the co-geodesic flow on the sphere may be written as 

1 



(6) 



R 2 sin 2 9 ' 



R 2 

P4> = 



cos 9 
R 2 sin 3 9 



The co-geodesic flow on the plane {z = 0} is equivalent to the above one through the stereographic 
projection. Since the metrics g and gs are in the same conformal class, the new system is easily 
derived using (j3]) with A = ( 1+c ° sfl ) _ 

However, on the plane we prefer to use the Cartesian coordinates (x, y) rather than the angular 
coordinates (<j),0). The latter arc related to the former by the transformation 



G 



arctan 



(f) 



9 = 2 arctan 



\/x 2 + y 2 
2R 



Hence, denoting by IR 2 the plane endowed with the metric gs, the Hamiltonian function for the 
co-geodesic flow on (R 2 ,gs) 

K geod : T*R 2 -> R : (q;p) ^ ^(G s p,p) 

is explicitly given by 

K geod (x,y;p x ,p y ) = l(x,y) [a(x,y)p 2 x + p 2 y ]. 

where q := (x, y) are the positions and p := (p x ,p y ) are the momenta. In order to derive this 
expression we set 

r i 2 

AR 2 + x 2 +y 2 

AR^/x 2 + y 2 

and we use the identities 



8R 2 



[AR 2 + x 2 + y 2 ) 2 ' 



smv = 



and 



ARy/x 2 + y 2 
AR 2 + x 2 + y 2 ' 

(1 + cos6») 2 _ 



sin0 = 



\A 2 + y 2 



8i? 2 



8i? 2 



(4i? 2 + a; 2 +y 2 ) 2 ' 

Note that a(x,y) corresponds to the term (sin#) -2 , while l(x,y) is just ^2^- To the Hamiltonian 
K geoc i is associated the Hamiltonian flow 



(7) 

where 

(al)(x,y) := 



I 



x = 2 (al)(x,y)p x 
is = 2l(x,y)p y 
p x = -{d x (al)p 2 x + (d x l)p 2 y ] 
k = -[d y (al)p 2 x + (dyl)p y \ 

32R 2 x 



2(x 2 + y 2 ) 



, 9 x l(x,y) = - 



{AR 2 +x 2 +y 2 Y 



, d y l(x,y) = - 



32R 2 y 



(AR 2 + x 2 + y 2 ) 3 ' 



3 Position of the problem 

We are now in the position to introduce a conservative force field on the sphere that perturbs, not 
necessarily by small amounts, the geodetic dynamics of a free particle governed by the equations 
discussed in the previous section. 

We place the singularity of the potential at the point Q = (0, R, R). As a naming convention, 
we shall often refer to it as the vortex point. On S \ {Q} we define the logarithmic potential U as 

U(P) :=-£log(||PQ||) 

where ||PQ|| is the three-dimensional Euclidean distance between P and Q, i.e. \\P Q\\ is the length 
of the chord between two points on the sphere. We denote by / the force field generated by the 
potential U, that is f(P) = VU(P). At any point P ^ Q it associates a force pointing towards Q if 
r > or in the opposite direction otherwise, and proportional to the inverse of the distance ||PQ||. 
Note that, unlike the planar case, the force field is not tangent to the manifold: at any point P one 
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could decompose the force vector into two components, one directed as the normal to the sphere 
the other tangent to the sphere. We take the first as balanced by the smooth constraint given by 
requiring that the motion happens on the spherical surface, therefore only the second contributes 
to the motion. 

The Hamiltonian augmented with the potential function 

H mech :T*(S\{Q})^R 

in the (<tj, (9)-coordinates is 

p 

H m ech{<l>,0,P<t>,P8) ■= H geod {4>,9,p rt>1 p g ) + — log(2i? 2 (l - sin 6> sin </>)). 

where 2i? 2 (l - sin0sin</3) = ||PQ|| 2 . 

On (M. 2 ,gs) the distance ||P<5|| 2 becomes the function 



b(x,y) := 

therefore we introduce the Hamiltonian 



2R 2 [x 2 + {y-2Rf 



4R 2 + x 2 + y 2 
K mech :T*(R 2 s \{0,V})^R 



r 

K mech (x,y 7 p x ,p y ) := K geod (x,y,p x ,p y ) + — log b(x,y) . 

We observe that logb € < ^ OC, (R 2 \{(0, 2R)}): indeed the point V = (0,2R) corresponds to the 
stereographic projection of the vortex Q £ S, while the origin is a singularity of the metric. 
Hamilton's equations associated to K mec h can be written as follows: 

x = 2 (al)(x,y)p x 
V = 2l(x,y)p v 

p x = -[d x (al)p 2 x + (d x l)p 2 + 17(8706(1, j/)- 1 d x b(x, y)] 



(8) 



Py = - [dy(a 1) P 2 X + (dyl) V \ + T / (8*)b( X , y)^ dyb(x, J,)] 

These equations govern the motion of a particle constrained on a spherical surface and subject 
to a force field generated by a logarithmic potential, as seen on a stereographic plane. 

Remark 3.1. Before we delve into the analysis of system (j5J) let us remark that when defining 
the potential function one could consider different notions of the distance between two points on a 
sphere. A reasonable choice could be the geodesic distance, that is the length of the shortest arc of 
a great circle passing through two points. More precisely, for any couple x,y £ S n (r) C the 
geodesic distance d§(x,y) is given by 

(9) d§{x, y) := r arccos — 5—, 



where < •,• > is the Euclidean scalar product in M. n+1 . In the following only the chord distance will 
be considered, but we guess that the local flow, that is the dynamics close and up to the singularity, 
should not be different when the geodesic distance is taken into account. On the other side, we 
expect the global flow to be slightly different. However a complete study of the dynamics with the 
geodesic distance, the differences and similarity with the dynamics on the plane and with the one 
here described, could be material for future investigations. 



8 



4 Energy hypersurfaces, regularization and flow 



We begin the analysis of system @ 
hypersurfaces associated to K mec h. 
by 



| with the description of the topology of the constant-energy 
For any fctl the hypersurface of constant energy h is given 



(10) 



{(%,y,Px,P y ) € T*X : K mech (x,y,p Xl p y 
(x,y,p x ,p y ) £ T*X : a{x,y)p 2 x 



h} 
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pi = 



l(x,y) 



h-r/(8n) \og(b(x,y)) 



where X := M 2 \{0, V} denotes the configuration space and T*X the phase space (the cotangent 
bundle over X). 

Since the function a(x,y) is strictly positive, for any value of h the motion is allowed only in 
those regions of the configuration space where the right hand side of the equation in the definition 
of T,h is positive (Figj2]). In the Lemma [4.11 below, the analysis is performed for T > 0: changing 
the sign of T simply switches the allowed region with the forbidden region. 




Figure 2: Zero-level curves of the function Eh(x,y) for different values of ft-. (T = R = 1) 



Lemma 4.1. For any fixed h let E h : M 2 — > R be defined by: 

(11) E h (x, y) := (h - r/(87r) log(6(x, y))) , 

and let hi = ^ log(2i? 2 ) 7 h 2 = ^ log(2i?). Then 

1. for every h > h 2 the surface Eh(x,y) is positive for any {x,y) £ M 2 ; 

2. for any h G (h\, h 2 ) there exists a disk D\ in the {y < 0} half-plane and containing the point 
(0, —2R) such that Eh(x, y) is positive for each (x, y) e M 2 \Dj l and negative otherwise; 

3. for any h < h\ there exists a disk Z? 2 in the {y > 0} half-plane, containing the point (0,2R), 
such that Eh(x,y) is positive for all (x,y) G D\ and negative otherwise. 

Proof. The logarithm is a monotone function, therefore the topology of the level sets of Eh{x,y) 
only depends on the level sets of b(x,y). For any 5 the 8- level set of b(x,y) is given by the point 
(x, y) lying on the circle Cs with equation 

The center of Cs is placed in the point Og = (0, 2 i$-& ) anc ^ the radius is r$ = r^JrLg] V^8R 2 — 8 2 . 
It follows that the function b(x,y) only admits values 8 in the range 8 £ [0,4i? 2 ]. For 8 = 4i? 2 
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the 5-level set restricts to the point P — (0,— 2R) while, as 6 decrease towards 2i? 2 , the level set 
consists of a circle completely contained in the {y < 0} half-plane. Moreover it can be checked that 
\Os — (0, — 2R)\ < rg, meaning that the point P is always surrounded by these circles. The value 
8 = 2R 2 is a singularity for the topology of the level sets: indeed the center of Cs as well as the 
radius r$ diverge. Note that b(x,0) = 2R 2 for any x and \im x 2 +y 2_> oc b(x,y) — 2R 2 . Then as 8 
decreases below 2R 2 towards zero, the circles C$ live in the positive-?/ halfplane and shrinks around 
the point Q = (0,2R). 

In order to develop a McGehee type transformation, we define the functions 



□ 



(12) 



(pi(r) := r e 1 l r 
¥>2(r) ■= l/r 



and, following the notation of jStFo03j , we introduce the change of variables 
(13) 



x = tpx(r) si 
y = yx{r)s 2 + 2R ' 



Px = y>2 (r) z x 

Py = ¥2{r) Z y 



where s = (si,S2) = (cos a, sin a) £ S 1 is a point on the unit circle. It readily follows that 

2 



a(r, s) 



l(r,s) 



8R 2 + 4R(p 1 (r)s 2 +<pi 
AR^/yi + AR? +4R(p 1 (r)s 2 

8R 2 



b(r, s) = 



2R 2 Vl {rf 



and (a i)(r, s) — 



8R 2 + <p 1 (r) 2 +4R<p 1 (r) s 2 ' 
1 



1R 2 + V 2 + 4R^ 1S2 ) 21 ' ■ v a.v 2{AR 2 + V 2 +AR Vl s 2 y 
hence in these new coordinates the energy surfaces can be written as 

2 

(r, s, z x , z y ) el+xS'xM 2 : a(r, s) z 2 x + z 2 y = j^—^ (h - r/(8vr) log(6(r, s)) 

We also observe that 

• lim r _j. + a ( r i s ) = 1 uniformly with respect to s; 

• lim r _>o+ b(r, s 



uniformly with respect to s; 
• lim r ^ + l{ r i s ) = l/(8i? 2 ) uniformly with respect to s. 



By taking into account the definition of the functions (fj, the right hand side of the equation 
defining the level set E/, reduces to 



E(h,r,s) := 



l(r,s) 



h - r/(87r) log(2i? : 



22-2/r' 



)+r/(87r)log(c(r,a)) 



where c(r, s) := 8i? 2 + ^i(r) 2 + 4i?Vi(r) s 2 = 8i? 2 + r 2 e _2 / r2 + ARre^ 1 / r s 2 . We observe that 

2Ti? 2 



(14) 



lim E(h, r, s) 



thus, as already implicit in Lemma 14. 1[ in the attractive case (r > 0) the vortex point lies in 
the allowed region of every energy level h, while the opposite holds in the repelling case (r < 0). 
However, a first important consequence of the change of variable above introduced is that in the 
variables (r, s) the kinetic energy remains bounded when a collision occurs. 
From now on, we shall only consider the attractive case. Thus we assume 

T > 0. 

The intersection between one (and hence every) energy hypersurface E^ with r = is called 
total collision manifold. In virtue of the limit (|14[) . we may conclude that 
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(a) (b) (c) 

Figure 3: Graph of the function E and of E = in the three cases: (a) h > fi2, (b) hi < h < /12, 
(c) h < hi. 

• the total collision manifold does not depend on the fixed energy level h; otherwise stated it is 
a boundary of every energy level; 

• it is diffcomorphic to the two dimensional torus T := S 1 x S 1 . 

From the dynamical viewpoint an important role is played by the zero set of the function E: 

Z h :={(r,a) EX : E(h,r,a) = Q}. 

In the following we refer to this set as the zero velocity manifold in S^. Rephrasing the results of 
Lemma 14.11 in terms of the new coordinates (r, a) it readily follows that Zh is empty in the first 
case (h > /12) and non-empty otherwise. (Figure [3]). In the second case, (hi < h < h>i), the zero 
velocity manifold is represented by a simple closed curve homeomorphic to a circle (or to a point in 
the limit h — > h^). The motion is forbidden in the region bounded by the curve. In the third case, 
the zero set can be seen as the graph of a single-valued function a 1— > r(a) and the function E is 
positive for < r < r(a), which is the region on the left of the curve shown in figure (J4]d) 




(a) (b) 



Figure 4: Zero Velocity manifold (a) in the second case h € (hi, h 2 ), (b) in the third case: h <h\. 



Regularization and McGehee coordinates 

We now use the new variables r, a and z in the equations of motion ([5]), In order to preserve 
the continuity of the flow with respect to the initial data, we need to ensure that the transformed 
system has an everywhere differentiable vector field. To this purpose we rescale the time variable 
in terms of the distance from the singularity with the effect to exponentially decrease the velocities 
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near the singularity. As a result the collision solutions (which are singular in the old coordinates) 
move along smooth orbits that asymptotically converge to the collision manifold. 
Let us define dr = ip2(f)<fii (r) dt and use the notation 

(z,s(a)) a '■= a(x, y)z x cosa + z y sin a. 



With the help of the identities 

<Pi(r) _ r 3 9?i(r) 
(p'i(r) r 2 + 2' 

the Hamiltonian equations in ([5J become 



<pi(r)ip' 2 (r) = r 2 
<P2(r)ip[(r) r 2 + 2 



(15) 



dr 
~aW 

da 

7h 

dz x 

~dV 



dzy 
~d7 



2r 3 



l(r,a) (z, 8(a)) a 



(2 + r 2 ) 

= l(r, a) (z y cos a — a(r, a)z x sin a) 
= - re~ 1/r2 [(al) x (r, a)z 2 + l x (r, a) z. 



-l(r, a)(z,s(a)) a z x 



re 1 / r2 [{al) y (r,a)z 2 x +l y (r,a)z 2 } 
-j— -l(r, a)(z,s{a)) a z y 



Stt 



r 

8^ 



r 3 e -l/r 2 ^M + 



b(r, a) 



b(r, a) 



+ 



where the subscripts in (al) x (r, s), (al) y (r, s) (resp. b x (r, s), b y (r, s)) denote the partial derivative 
with respect to the old cartesian variables. The derivate function is then evaluated in the new 
coordinates at the point (r,a). The equations above are no longer singular at r = 0: in fact, by 

computing -j-(r, a) one finds -y-(r,a) ~ (pV 1 as r — > 0. Thus the time change produces the effect 
oo 

to regularize the singularity. In addition the {r = 0} manifold results to be invariant. 

From a naive point of view, the study of the flow on the collision manifold could appear mean- 
ingless, since the manifold is the image of just a singular point where the orbits cease to exists. In 
reality, the properties of the flow on such manifold yield informations on the behavior of the orbits 
close to the singularity. 

In order to simplify system (fT5")) , we introduce a further change of coordinates and time rescaling. 
Using the energy relation 



(16) 

let us define tp an d cr such that 



a(r, a) zl + z 2 y = E(h, r, a), 



(17) 



z x = y E(h,r,a)/a(r,a) cosi/j, 



z y = J E(h, r, a) sin^ 
dr = y E(h, r, a) da 



Let us denote with A±(r,a,z), A^ir, a, z) the right hand side of, respectively, the third and fourth 

equation in (fl5)) . so that — — = A\, — A 2 . Therefore, on any fixed-energy shell, the system 

dr dr 



12 



given in (|T5|) reads 

dr 2r 3 E(h, r, a) 



(18) 



I (r, a) ( y/a(r, a) cos ?/> cos a + sin ?/> sin a) 



da (2 + r 2 ) 

E(h,r,a)l (r, a) (sin ^ cos a — ^/ a(r, a) cos ip sin a) 

; S(r,a,V) 



fid 



d?A 



where -B is given by: 

B(r, a, tb) := —y/a sintbAi + y/asmip cos ip-^j- [ \ — I + A 2 cos?/; ^-(V E) sin ?/> cos ?/>. 

ar I V a / ar 



Flow and invariant manifolds 

Recalling that £ is everywhere positive, the restpoints of (fT5|) correspond to solutions of the following 
systems: 



(19) 



where 



r = 

f 2 (r,a,ijj) = 
B(r,a,if}) = 



fi(r,at,ip) = 
f 2 (r,a,tfj) = 
B(r,a,ip) = 



E(h,r,a) = 
B(r,a,?/>) = 



/i(r, a, -0) := cos?/; cos a + sin?/' sin a, /2(»™) a, ?/>) := sin?/; cos a — ya cos?/' sin a . 

We immediately discard the second system as it allows no solutions. We note that a — > 1 for r — > 0, 
thus the first system reduces to 

r = 

sin(?/> — a) = 
B(r,a,ip) = 

whose solutions correspond to fixed points on the collision manifold. The existence of solutions of 
the last system depends on the energy level h: if h > h 2 the zero set of E is empty and no solutions 
exist. For h < h 2 , some solution may exist. 

Summarizing, any restpoint either lies on the collision manifold or on the zero velocity manifold. 
Let us first consider the collision manifold: the asymptotic analysis of the function B(r,a,ip) on 
the collision manifold (see Appendix A) gives 



lim B(r, a, ip) = 0. 

!-->0+ 



It follows that 



Lemma 4.2. The equilibria of the vector field given in (1181) lying on the total collision manifold 
consists of two curves. In local coordinates (r, a, ?/;) these curves are given by 

(i) &> x = (0,a,a); 

(ii) ^ 2 = (0,a,7r + a). 

Proposition 4.3. For each a, the equilibrium points 

(0, a, a) £ 

and 



(0, a, 7T + a) £ @> % 



are degenerate saddles. 
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At ^ 2 
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Table 1: Dimensions of the invariant manifolds along the equilibrium curves ^\ and 

1. dxmW u {^i) = 1, dimW a (^i) = 1, dimW ^) = 1. 

& dimW u (^2) = 1, dimW 8 (^2) = 1, dimFy°(^ 2 ) = 1. 
Proof. The flow on the collision manifold is given by 



(20) 



da T \ 
— = — sm(V> - a) 
da 47r 



dip 
da 



= 0. 



whose orbits are parallel to a-axis and flow from to 8?\. The stability of the restpoints is 
determined by the eigenvalues of the Jacobian matrix of (IT8l) . It follows (see appendix [2| that for 
any point P\ £ and Pi € ^ 2 the eigenvalues are 

x r = o r A r = 

Pi e &>i { X a = ~ , P 2 G ^ 2 =>• J A Q = -J- 

47T 47T 

A^=0. [ A v , = 0. 




3fr/4 



Figure 5: The collision manifold, the curves of restpoints , 
of equations ([2"0|). 



(blue) and ^ 2 (red), and the vectorficld 



and are coherent with the dynamics restricted on the collision manifold as given by (|20[) . where 
&i is an attractor and is a repeller (figure [5]). However, the presence of null eigenvalues implies 
that the linear approximation of the flow, taken alone, does not provide enough information to 
determine the qualitative dynamics close to the equilibrium points. As flow in the tp direction is 
null (and in fact ip can be regarded as a parameter for an equilibrium point), in order to determine 
the asymptotic behavior close to P\ and P2 , it is enough investigate the dynamics restricted to the 
(r, a)-plane. The proof that the two equilibrium curves are indeed degenerate saddles follows by 
direct integration of the system once the equations have been expanded around the equilibrium 
point in Taylor series. We omit the details and we refer to the equivalent proof of Lemma 7.4 in 
|StFo03j . □ 
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Definition 4.4. We shall say that the flow on the collision manifold is totally degenerate if the 
unstable manifold of an equilibrium point Pi £ £P\ , coincides with the stable manifold of some 
equilibrium point P 2 £ 2P 2 . 

Lemma 4.5. The flow on the total collision manifold is totally degenerate. More precisely 
(1) W U (P 1 ) = W S {P 2 ); 
(11) W U (P 2 ) = W s (Pi); 

where Pi € S?\ and P 2 G S? 2 are chosen in such a way the last coordinate of the two points agrees. 

Proof. The proof of this result follows by a straightforward integration of the equations of motion 
on the total collision manifold r = 0. □ 
A direct consequence of the previous result is the following: 

Corollary 4.6. (Existence of heteroclinic connections,) There exists an heteroclinic connec- 
tion between each equilibrium point P\ € S?\ and the point P 2 G 3P 2 where P\ , P 2 were chosen in 
such a way that they have the same projection on the first and third coordinate. 

Proof. The proof of this result follows immediately by the previous result. By the fact that r = 
and ip is constant, it follows that the non equilibrium solutions are in the (a, ip)-p\a,iie lines parallel 
to the a-axis. Moreover each point of equilibrium on &\ is attracting while each equilibrium point 
on S? 2 is repelling. □ 



Moving out of the collision manifold, the two lines &\ and ^2 exhibit the opposite stability 

dr dr 
character: indeed — > for ip = a and — < when ip = a + w, meaning that the system goes 

da da 
into the collision along £? 2 and escape from the collision along 8?\ . 

Next we examine the restpoints and the flow on the zero velocity manifold. This is more easily 

accomplished by looking at the system given in (|15p . Restpoints, in fact, are not changed by a time 

scaling. Since the zero velocity manifold coincides with the zero set of the function E and by taking 

into account 

Setting E = in the energy relation (1161) . which implies z x = z y = 0, it follows that on the zero 
velocity manifold the dynamical system (1151) reduces to: 

dr 

d~r=° 



(21) 



da 

dz x 
dr 

dZy 
1^ 



_^ r 3 e _ 1/r 2 b x {r, a) 



8tt 



b(r, a) 



r r 3 c -i/r 2 b y( r ^) 
8n b(r,a) 



The restpoints on the zero velocity manifold (if any) correspond to the solutions of the equations: 

b x (r,a) = by{r,a) = Q 
b{r, a) b(r, a) 

An elementary calculation shows the following result. 

r 

Lemma 4.7. For h — h 2 (defined above as h 2 :— — log(2i?)J there exists only one restpoint on the 

An 

zero velocity manifolds at P := (r*,37r/2) for <^i(r*) = 4i?. For h 7^ h 2 there are no restpoints. 
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5 Global flow and dynamics on the sphere 



It is now possible to bring back on the sphere the results found on the stereographic plane in the 
previous sections. 

In terms of the coordinates (6, 9) on the sphere the Hamiltonian reads as: 

H(9,6,p e , P<f ,) = ^ (^~e P * +p e) + £ lo s( 2i?2 ( 1 " sin sin 0)). 

Note that the vortex is located at (<j>,0) = (tt/2,tt/2), therefore for ((f), 9) — > (tt/2,tt/2) the dynam- 
ical behavior becomes unknown since the vectorfield ceases to exist. 

Let us call vortex half-sphere the half-sphere centered around the vortex point, and antivor- 
tex half-sphere the complementary half-sphere; let us call vortex-parallel any circle on the sphere 
equidistant from the vortex and vortex-meridian any great circle passing through the vortex. Fi- 
nally let us call antipodal point the point on the sphere opposite to the vortex point. On the sphere 
the results of Lemma T4. II can be rephrased as follows. 

Theorem 5.1. If h < hi the motion is allowed in the region of the sphere containing the vortex 
point and bounded by a vortex-parallel that lies on the vortex half-sphere for h < h\ and in the 
antivortex half-sphere otherwise. If h > hi then the motion is allowed everywhere on the sphere. 

Moreover lemma |4~T1 is rephrased as: 



Theorem 5.2. Fork 
point. 



hi the zero velocity manifold consists only of one point which is the antipodal 



In order to understand the global dynamics it is useful to show the existence of a second conserved 
quantity, analogous to the angular momentum for planar dynamics. To this aim it is convenient 
to move the vortex point at the north pole N — (0, 0, 2R) (or, equivalently, to redefine the pa- 
rameterization of the sphere). Obviously this does not change the dynamics. Note that the curves 
{9 = const} and {</> = const} now respectively correspond to the vortex parallels and to the vortex 
meridians. 

In this setting, the dynamical system reads as 

■ 1 



(22) 



R? sin 2 ( 



Pe_ 

R 2 

= 



s'm9 



Pe 



i? 2 sin 3 <r 87r2i? 2 (l +cos6») 



corresponding to the Hamiltonian 



H(9,6,p e ,p^ = ^ (-r^Pl + Pe) + £ log(2i? 2 (l + costf)). 



Going to the Lagrangian formulation, we can write the Euler-Lagrange equations 

^(R 2 j>sm 2 9) = 
dt 



R 2 9 - R 2 sin cos 6^ 2 



Btt 2i? 2 (l + cos( 



0. 



The first shows the existence of a conserved quantity, namely the spherical angular-momentum 
/ = R 2 sin 2 66. It follows that 
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Figure 6: Particular orbits lying on the vortex-paralles and vortex-meridians 



Lemma 5.3. A necessary condition for a solution to either collide with the vortex or to reach the 
antipodal point is I = 0. 

Proof. Writing the energy relation in terms of (<f>,9,4>,6) and substituting <f> = R2 . 5 e , it follows 
that a solution exists only for those 9 satisfying 

2R 2 hsm 2 8- —R 2 sin 2 6>log(2i? 2 (l + cos 6)) - I 2 > 

4"7T 

Recalling that the vortex is placed at 6 = tt, and that the antipodal point is at 6 = 0, it follows 
that if either a collision occurs, or the point goes to the antipodal point, then I 2 < 0. □ 

Looking at the system (f2"2")l , one can easily prove the existence of particular solutions as depicted 
in Fig (5] 

Lemma 5.4. (i) Any vortex-parallel on the vortex half-sphere is the support of a periodic orbit. 

(ii) The vortex-meridians are flow-invariant. 
Proof. For any 6 £ (§,7r) the curve 

78 (t) := (cf>,e,p^p g )(t) = (fo+t jpff 0)a » g 'P«'°) 

with p^ 2 = — 2 ctJa"^^^ g) is a solution of the system. Note that the previous relation can not be 
satisfied if G (0, 7r/2], which implies that only the vortex-parallel placed in the vortex half-sphere 
are support of periodic orbits. Moreover the period of 70 (t) tends to zero as 6 goes to ir/2 or 
7T. This proves statement (i). Statement (ii) immediately follows by noting that any initial data 
(4>,9,p^,pg)(0) = ((jP ,9° ,0,pg) leads to an orbit traveling on the {(f> — 0o} vortex meridian. □ 

A consequence of the previous lemmas is that any orbit with energy h > hi that passes through 
the antipodal point will end into the collision. 

The existence of heteroclinic connection anywhere on the total collision manifold of the regular- 
ized flow fCorollarv l4.6[) provides the way to extend beyond the collision the orbits of the singular 
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flow. In fact, if ^(t) = (4>,8)(t) : [0,T S ) — >• S is a collision solution ending in the singularity at time 
T s , we define the collision-transmission solution as the path 7 : [0, 2T S ] —> S as 



where (0y , 9y) are the coordinates of the vortex point. The extended flow obtained by replacing the 
singular trajectories with the collision-transmission solution results to be continuous with respect 
to the initial data. The same result for a single logarithmic center on the plane has already been 
proved, among others, in [CaTej with a completely different technique. 

Finally, from the above discussion it follows that the collision-transmission solution behaves 
in three different ways, depending on the energy level h: if h < I12, after the ejection from the 
singularity, the particle reaches the zero velocity manifold, then it reverses the motion and falls back 
into the vortex point; if h = /12 after the ejection the orbit asymptotically reaches the antipodal 
restpoint (this is an heteroclinic orbit between a point of and the single restpoint on the zero 
velocity manifold), and if h > h,2 after the ejection, the orbits travels along a vortex meridian, 
passes through the antipodal restpoint, continues the motion on the opposite meridian and falls 
down again into the singularity. 

A Useful asymptotics 

In this appendix we list some asymptotic limits of the functions appearing in the equation of motion 
in McGeheee coordinates. They are useful to compute the spectrum of the eigenvalues associated 
to the fixed points. 

On the total collision manifold r = 

All the limits below are computed with respect to r and are uniform with respect to the other 
variables. For the function E, we have: 




lim E(h, r, a) 



AT i? 4 



lim d r E(h, r, a) = lim d a E(h, r, a) = 



7T 




For the functions a, b, we have 



lim a(r,a) = 1, lim d r a(r,a) — 0, lim d a a(r,a) = 

r->0+ r->0+ r->0+ 




lim b(r, a) = 0, lim d r b(r, a) = 



For the functions a x ,b Xl a y ,b y , we have 




lim a x (r,a) — lim a y (r, a) — lim d r a x (r, a) = 




lim b x (r, a) = lim b y (r,a) = lim d r b x (r, a) = 





lim d a a x (r,a) —0 lim d a a y (r, a) 








lim d a b x (r, a) — lim d a b y {r,a) 









lim 



r 3 e 1/r2 b x (r,a)d a b x (r,a) 



r 3 e 1/r2 b y (r,a)d a b y (r,a) 
b 2 (r, a) 



= 



6 2 (r, a) r^o+ 



IS 



For the functions z x , z y , we have: 

/r [r 
— 2R 2 costp, lim z v (r, a, ip) = \ — 2R 2 simp, 
7T r^0+ V 7T 

lim d r z x (r, a, ip) = 0, lim d a z x (r, a, ip) = 0, 

■r^0+ r->0+ 

/r /r" 
— 2_R 2 sin^, lim d^z y (r, a, tp) = \ — 2R 2 costp, 
7T l — >0+ V 7T 

lim d r z y {r, a, tp) = 0, lim d a z y (r, a, ip) = 0. 

r— f0+ r— >0+ 

As consequence of the above asymptotic behavior it follows that 

lim Ai(r,a,ip) = 0, lim A2(r, a, ip) = 0, lim £?(r, a, ip) — 

r->0+ r->0+ r->0+ 

lim d a Ai(r, a, ip) = 0, lim Sq,^^, a, VO = 0. 
Denoting by J := (Jij)ij the variational matrix on the total collision manifold it follows that 
Jn = 0, Ji2 = 0, Ji3 = 0, J 32 = 0, J 33 = 0. 

r 

J21 = 0, J 22 = -— cos(V> - a), J 23 = -J22, 

lim c^ T a(r, a) = 0, lim d 2 aT E{h, r, a) = 0. 

The limit involved in the computation of the term J31, in general may not exists. However at the 
restpoints ip = a or ip — a + tt this limit actually exists and this implies that J31 = 0. 
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